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Abstract 

Accurate approximations are presented for the time development of 
both edge conditions and internal structures of a blast wave with 
shock heated electrons, and equal ion and electron temperatures at the 
shock. The cases considered evolve in cavities with power law ambient 
densities (including the uniform ambient density case) and have 
negligible external pressure. Account is taken of possible saturation 
of the thermal conduction flux. The structures evolve smoothly from 
those given by Cox and Edgar (1981) to the adiabatic structures given 
by Cox and Franco (1981). 
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This work continues a project begun in Cox and Edgar (1983, 
hereafter Paper I), modelling spherically symmetric supernova blast 
waves in a power-law ambient density distribution. The gas is assumed 
to be ideal, totally ionized, and the helium abundance is 10% of the 
hydrogen abundance by number. We have assumed also that non-coulomb 
processes heat the electrons in the shock, so that the electron and 
ion temperatures are equal immediately behind the shock. Thermal 
conduction is then important at early times when the temperature is 
high, but less important later, when the electron temperature and 
conductivity are lower. Our formulation of the thermal conduction 
flux follows that of Cowie (1977), which accounts for saturation 
effects. 

We report here an approximate solution to the evolution of 
temperature, density, and pressure structures of such a supernova 
remnant. We separate, via physically reasonable approximations, the 
immediate postshock behavior from the evolution of the system as a 
whole, and solve for the evolution of the postshock conditions.' We 
then separate, again by reasonable approximations, the dynamical and 
thermal evolution of individual gas parcels from that of the overall 
structure, and follow conditions in these parcels. Knowledge of the 
history of each parcel would facilitate the calculation of ionization 
fractions, which may be far from equilibrium. 

We thus continue a study of the first order effects of thermal 
conduction on . the structure, spectrum, and surface brightness 
distribution of a ”non-radlative" remnant, including as much physics 
as possible, while avoiding the numerical solution of coupled partial 
differential equations. 
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The notation for the variables we have used is listed in Table 1, 
along with that of their post-shock derivatives ( 3£n f/ 3£,n R) and 
normalized forms. As in Paper I, the jump conditions and their 
derivatives are given by 
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p* = x*-a). 


( 2 ) 


where p « was assumed for the ambient density, the asterisk 

represents the logarithmic deriva-tive f* H 3£n fg/ 3£n Rg, Fg is the 
thermal conduction flux just inside the shock, and v^ is the shock 
velocity. In addition, ambient pressure has been neglected, and 
Pg = -3 has been assumed. This last is equivalent to assuming PgRg is 



-A- 


constant, which cannot be exactly correct, but reflects the general 

effect of conservation of energy. In Paper 1 it was shown that the 

early and late asymptotes have p„rJ which differ by only a few percent 

s s 

at most (depending on m). This assumption may tend to smooth over 
real short term excursions in pressure but prevents spurious 
approximation-induced variations and thus provides a physically based 
stability to the behavior. 

Writing the post-shock equations of motion for a gas parcel in 
terms of the logarithmic derivatives given in Table 1 produces several 
relations among these quantities. Mass and momentum conservation give 


e = 


^s-1 


^s^s 

^s-1 
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= (B - l-a) + Xg(3 + 


3 ^ F, 


The energy equations for the ions and electrons can be separated 
(Paper I) and written in the form 


Dt '' n ^ ^ n p 


(4) 


ff— -g)Tp/2 

Dt^ n ^ k.Z-'k.Z 


( 5 ) 


where T = (I.IT- + 1.2T )/2.3 is the average temperature, g = T /T, 
and X. = (1/r2)(3/3R)(r2f) is the divergence of the thermal conduction 
We have assumed here that T^ » Tj|^/20CK>i 


flux. 


The cons t ant 
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a ^ In A/153 cgs where In A = iln[1.2 x 10 5 . Also, 

n = Hy + rijjg = I.Idjj, P£ = nkTj^, = (1 .2/1 . l)nkTg , so that the ideal 

gas law becomes p = (2.3/l.l)nkT. 

Equation ,(5) can be solved at R = R (using the assumed g„=l 

S' s 

appropriate to electron and ion temperature equilibration by 
non-coulomb processes in the shock) for v = (9£.n g/ 9f,n R)^ to produce 


1.1 ^"^s f^S^S-v 

" ■ TT -3- 

s 


( 6 ) 


Writing equation (4) in terms of n = ~ ( BZnT^/ dZn R)g and making 
use of the fact that -q=\H-(S-a), it follows that 


(5x« - 8)(n + v) = (6x| - llx„ - 4) - tJXg(4Xg - 7) + x„x*(5xg - 7) 




-(X, - 2)(4 - xj 


Rs^s 

S-' V * 
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which, together with equation (6) gives us q as a function of x^, x*, 
and Rglg/Fg. 

To evaluate the thermal conduction flux, we use the form given by 
Cowie (1977) for the ratio of true flux to classical (i.e. 
non-saturated) flux: 



4.6 nX 


( 8 ) 


- 6 - 


where (|)g is a plasma parameter of order unity (we have used ({ig = 1 
throughout), and X is the mean free path for electron energy exchange. 

In our case, for conditions just inside the shock front, the 
second . term in the denominator can be written (using 
X=m^^^ bT2/(1.31ngk^^2) where b = 6xl0~^ (30/iinA) cgs, from Cowie and 
McKee (1977), the jump conditions relating F^, Vg, Pg> and the 

ideal gas law) 

4.6 riX 1 ^“^s ^c .Q 


whence it follows that 

-Li = 1 - _J_ 

^c'’® 10<^s (Xg-l)l/2‘ 


( 10 ) 


Notice that the degree of saturation of the thermal conduction 
flux just behind the shock is simply related to x . The conduction is 
never heavily saturated and becomes classical as Xg approaches 4. The 
classical flux at Rg can be written 

3Tg 

IF 


Fc)^ - -b I|'2 




bl|'2 


( 11 ) 


Hence 
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T1 = 

bl7/2[i 
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In this equation, T_ and F„ can be eliminated using equation (1) and 

S d 

the ideal gas law. The result can be written 


x|(4-Xs) 

" Knotty- 


2 7-301 
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where z is a normalized radius defined by 


27-30) = t'17- 
'■R. ^ 
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The z^~^“ dependence follows from n^ « and p_ «c R~^. In this 
formula, m = p/n = (1.4/l.l)m^. The formula for determining the 
costant Rj^ depends on energy integrals which we shall discuss later. 

From equation (13), it is clear that n depends very sensitively 
on Rg in the vicinity of Rj^. For small values of z, q will be 
extremely small and the electron temperature will be flat. For large 
values of z, however, x^ is driven to 4 to keep n from becoming very 
large. In fact, q approaches 4-3 uj in the adiabatic limit. 

In order to reduce the evaluation of the dynamics of this blast 
wave to the solution of ordinary differential equations, one further 
piece of information is required, namely the behavior of the parameter 
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RsAg/Fg* This parameter was discussed at some length in Paper I for 
the early time asymptote. In particular we can write 

— = fr^Jidr = / r^xf— )dr = <-?> , (15) 

Rg i 0 'x' X 3-0) 

where <£/x> is the mass average of the energy loss rate per n^ 
particles due to conduction. If all electrons lose energy at 

approximately equal rates via this mechanism, then <il/x> ~ 

£gRg/Fs = Xg(3-o)). It was found in Paper I that in fact £/x varies 
very little over the structure of the young blast wave (see figures 
lb, 2b, and 3b of Paper I) and that RgJ^g/Fg = Xg(2.53-o)) fit the edge 

results very well. This is certainly the value we need at the 

beginning of our evolution, but we are' about to argue that it is 
probably an adequate approximation throughout. 

At late times (z just slightly greater than 1 constitutes late 
times) X is driven to 4 as we have seen. Since R_ £„/F„ is everywhere 

in the equation multiplied by (4-x ), the particular form of 

is unimportant. Thermal conduction simply fades away. 

During the transition period, however, there is a very complex 
set of interactions. Thermal conduction is suddenly finite, allowing 
significant gradients in T^ to develop behind the shock. At the same 
time (e.g. Cox and Anderson, 1982, hereafter CA) , coulomb heating is 
becoming important near the edge. The coulomb heating is proportional 
to T^-Tg, however. If thermal conduction is truly Insignificant, Tg 
will remain at Tj^^ behind the shock and coulomb heating will do 
nothing. If thermal conduction is still important, T^ will fall just 
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far enough below that coulomb heating can supply the losses. Thus 
a competition develops between two processes for control of the 
gradient in T^. At some moment, control will pass decisively to 
coulomb heating, and the difference between and will become 
imperceptible over a progressively larger region. This kind of 
behavior is shown in both CA and Cox and Franco (1981, CF) for blast 
waves with coulomb heating only. 

Because of the sensitivity of both coulomb heating and thermal 

conduction to the post shock gradient of T^ as one enters the 

transition, the conduction flux and T^^-Tg will be difficult to 

approximate reliably from first principles. The ratio i_/F„, however, 

s s 

can be expected to vary rather smoothly. Apart from a very short 
lived transient during which near equilibration has been achieved only 
at the immediate edge, the ratio should never be very far from the 
mass average. (The mass is far too heavily concentrated near the edge 
for the result to be otherwise.) In short, approximating 
throughout by Xg(2.53-uj) may provide a very slight smoothing of the 
behavior near the onset of the rapid transititon, but as a rule it 
should be very nearly correct; it makes no a priori judgement about 
the magnituides of F or it is exact at early times when thermal 

conduction is most important; and it does not interfere with the 
approach to negligibility of thermal conduction at late times. 

With this approximation, equations (6) and (7) combine to produce 
an ordinary differential equation x*(z,Xg); 
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dx. 


dz 


5TC5^ U 5 .,- 8 ) [n + ^ ^ X 3 C 2 . 53 -«) ] (16) 

s s 

- (6x|-llXg-4) + ajXg(4xg-7) + (Xg-2)(4-Xg)Xg(2.53-aj) 


with T| given in (13). 

The Initial conditions for this equation are provided by Paper I. 
In the early time (small z) limit, x*=0, and the n term can be 
neglected (thus reducing equation (16) above to equation (17) of Paper 
I), so that the solutions for the compression factor x^( (u) given in 
Table 1 of Paper I can be adopted as Initial values. 

It should be noted here that the Paper I results for x^( w) and 
that (^s^s^^s^ “ (2.53-o))Xg follow from the assumption that the 
interior density structure is approximately that given by Kahn (1975) 
as generalized by CA. This approximation (which automatically 
conserves mass) was tested against analytic results for adiabatic 
blast waves by CF and found to be accurate to within a few percent 
(and much closer near the edge where the mass is concentrated). 

We have found solutions to equation (16) by numerical integration 
for u) = 0, -2, and -4 (uniform density and cavities of various 
steepness). The integration provides x^ and x* as functions of z, and 
introducing these into the post-shock derivative equations supplies a, 
3, V, and n all as functions of z. When the results are combined with 
the jump conditions and ideal gas law, they completely specify the 
immediate post-shock values of p^, Ug, T, and Tg and their spatial 
derivatives. The results are shown in Figure 1. A good approximate 
method for generating these functions for application to a particular 
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problem is presented in section III. A few properties of these 
solutions are given in Table 2. 

Shock radius as a function of time follows from the first jump 
condition in equation (1): 


-1 


Pb = 


Po^i 




( 17 ) 


where e(m) = 2EQ/3VgPg is the energy integral given in Paper I, and 
repeated in Table 2 for reference. Here is the volume of the 

remnant. Defining a natural time unit of 


2ir e(oj)pQ(R )r5 

h ■ [ 1 


(18) 


we obtain the ordinary differential equation 

ll- (19) 

dz ^ Xs 

for the age of the remnant x B t/tj^ as a function of shock radius 
z = Numerical integration is straightforward. The results are 

well approximated by a power law x « the maximum discrepancy 

amounting to about 3%. Figure 2 shows the behavior predicted by this 
integration. 
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II. The Interior Structure 


In this section we will derive the approximate internal density, 
pressure, and temperature structures of the remnant. The approach we 
will follow is to solve the equations of motion for selected gas 
parcels from the time they are shocked as they move into the interior 
of the remnant. 

A convenient coordinate to use in tracking gas parcels is the 
mass fraction y = M(R)/M, where M(R) is the mass enclosed within a 
radius R (which is constant for a given gas parcel, though R changes) 
and M is the total mass of the remnant, which Increases as the shock 
encompasses more material. 

The density structure will be assumed (as above) to be given by 
the approximation of Kahn (1975) as generalized by CA: 


X = 



- l)r‘ll 


y(r) = {r5/2 exp[[xg - . 1 ]}^ “ 

Z q 


Xg(3+a) - (3-0))x| 



( 20 ) 


The approximation conserves mass, provides the correct density and 
slope at Rg, and the correct logarithmic slope at R<<Rg. 

For a given gas parcel 1 which was shocked at we can 

calculate its present radius r = R/Rg by noting that yj^ = (zj^/z)^~“ 
(since M(R) has increased in._ proportion to z^~“), and solving 4:he 
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second of equation (20) numerically for r. Using (20) we then have the 
local density, since n= x n^(z=l) z~^. 

To find the remaining quantities of interest (p, T, and T^), we 
need two equations (plus the ideal gas law). These equations are 
provided by equations (4) and (5) above. It is convenient to choose 
as variables the "adiabatic constant” K, conveniently normalized, and 
another quantity C: 


K 


1 t3/2 

A n ’ 


- ( 2-3 1 

= [l3- g] 


3/2 


K. 


A 


1.4 ('1.4>,3/8 
1.1 ^ 2 . 3 -* 


, Ps<2=l) 

^rPo(z=l)]!>^3 


} 


3/2 

9 


( 21 ) 


where m^ is the mass of a hydrogen atom and g=Tg/T as before. The 
normalization constant A is chosen so that the initial value of K is 
particularly simple in form; see below. Then using the first 2 
members of equation (1), equation (6), and replacing time derivatives 
with z derivatives using z(D/Dz) = (R^/Vg) (D/Dt) we can rewrite 
equations (4) and (5) in the more convenient form 


DK ^ _ Kr£ X 
Dz z 2 . xj 


12 

11 


JXM], 


( 22 ) 


and 
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DC 

Dz 


-Gflll 

'^ 1.2 


- g) 


3-(u 

1/2 (1-g) „ 2 




(23) 


where y = p/p^ = (x/Xg)5/3(K/Kg)2/3, g = ± (^)tr H-37, and 
x(y) = ( ^/p)/( ^e/ Pc) i-S the normalized divergence of the thermal 
conduction flux. 

Using equation (16) from Paper I and the approximation discussed 
above that Rg^.g/Fg ~ Xg(2.53-o)), a bit of algebra produces the result 


a 1-2.3 


" 2.3(2.'53--^)(-4-xT 


(24) 


where x' = 3x/ 8r and r is given in terms of y by inverting equations 
(20). This formula was derived assuming self-similarity, and is 
therefore only valid at early times (the situation addressed in 
Paper I). We will take the same functional form for x(p) ^t all 
times, which necessitates using early-time asymptotic values for x 

s 

and q in equation (24) and the supporting equations (20). Thus at 
early times when the thermal conductivity is large it is correctly 
accounted for, and while the approximation is less accurate at later 
times, the right-hand side of equation (22) is approaching zero since 
v«(4-Xg). Thus the approach to the final adiabatic (K«<r372/jj 
= constant) condition is preserved, though the rate of approach may 
not be correct. 

We now have a pair of coupled equations, which can be solved - 
subject to the initial conditions (i.e. post-shock values) derived 
above, rewritten in the form 
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Ks = ^ 


5/2 2 ^ 


C,- [^- 


( 25 ) 


The solutions to these equations for representative parcels has 
been carried out, and snapshots of the resulting structure are 
presented in Figures 3, 4, and 5. The density plots are normalized to 
p^(z=l) or nQ(z=l), the pre-shock density at Rg=R^. The temperature 
and pressure graphs are plotted in units of the post-shock values at 


51 


. 1.03x10’ K 


1 ^12.2 pc 


Xg(z=l)nQ(z=l) R^ 


(26) 


Pg(z=l) = 2.97x10 9 dyne cm ^ , 


where e((o) is the energy integral from table 2, is the explosion 
energy in units of 10 ergs, and R^ can be calculated from 
equation (27) below. 
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The reader will notice the small waves near the edges of the last 
snapshots in the temperature and pressure plots given in Figures 3-5. 
We expect the sudden transition to introduce such waves, but the 
particular form shown is probably rather heavily influenced by our 
approximation scheme, in particular the constraints of the density 
structure to the Kahnian form and the post shock pressure to 
development. Considering that the model was probably forced to evolve 
somewhat more smoothly than the actual transition, the wave amplitudes 
should probably be regarded as a lower limit to both the uncertainty 
on local conditions and the expected transient amplitudes. 


III. Application of These Results 

In general one wishes to calculate the conditions in a remnant 
with some particular explosion energy, E , shock radius R„, and 
pre-shock density Pg(Rg)«R~“ and one is hopefully content with one of 
the values u = 0, -2, or -4. If one wishes to assume that there is no 
non-coulomb heating in the shock or no thermal conduction, one turns 
directly to CF . If one wishes to include a significant external 
pressure, is content with oj=0, and either no non-coulomb heating or no 
thermal conduction, one turns to CA. Those two papers are also useful 
if one is interested in an epoch significantly later than the electron 
equilibration time (discussed in either paper). If one is interested 
in a case with non-coulomb heating (T^=Tj^ at the shock) and strong 
thermal conduction for times much less than the equilibration time, 
one turns to Paper I. Only if one is forced to be interested in the 
structure of a remnant with non-coulomb heating and thermal conduction 
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in the vicinity of the equilibration epoch will one wish to follow the 
procedure below. 

For given values of E^, p^CRg) and o), one first finds 

e(o))=2EQ/3VgPg from Table 2. One can choose a value for the early 
asymptote, the late asymptote, or somewhere in between, depending on 
whether the time of Interest is relatively early, late, or 
intermediate in the transition. (The values are very similar and one 
should worry very little over the choice.) From this choice, one can 
then calculate p^. The mean mass per nucleus is m=1.4mjj/l.l and 
no=Po/m- 

The radius and time units and t can then be calculated from 


R = 12.2 pc 

1 ^ ^^ 12.2 PC'* '■eCw)'' 


■ e( (u) ■' £,nA 


[no(Rs)P> 


-3 i7-3o) 


(27) 


t = 3.19xl0V {(- 


E 3+u) 
51 


^ 5- CO 1 

^-4(0 , 51 ~ 2 - ,30 r , . i-4 .7^ 

12.2 pc^ ^eCco)-* '■ilnA'* 1- J ^ 


which follow from equations (14), (17), (18), the ideal gas law, and 
the fact that n (R The present maturity of the shock is then 
Zp=Rg/R^, and the age of the shock can be read from Figure 2. The 
postshock temperature and pressure and the preshock density at z = 1 
(which are the units of the ordinates in Figures 3-5) can then be 


found from 
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PoCRi) 


PoC^s) 



(28) 


and equations (26). 

If one is interested in only a cursory picture of the remnant state, 
the postshock compression can be read from the appropriate figure. 
The shock velocity, postshock mass velocity, and thermal conduction 
flux are found from equation (1) and the postshock temperature (and 
electron temperature since = T was assumed at the shock) from 
T = 1.1 p/(2.3 nk) . The logarithmic derivatives of p, p, T, T^, and 

T /T, can all be read from the graphs. If z is close to one of the 

e p 

those represented in Figure 3, the remainder of the normalized 

structure is shown. 

If, however, one requires a rigorous numerical description in 

order, for example, to calculate the past evolution to and the 

ionic concentrations now present, this is most easily carried out by 

using an analytical approximation to our calculated behavior of x . 

s 

This approximation must be exceedingly accurate to allow the 

derivation of the other parameters. We have found that writing x^ in 
the form 


1 + ciW + coW^ + CoW^ 

Xs = Xg(0) { _} (29) 

1 + dj^W + d2W^ d3W^ 


where the coefficients are as given in Table 3 and W = z^ provides 
sufficient accuracy. Then x* = (7-3(o)W(dXg/dW) /Xg is straight- 
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forward. The postshock, quantities then follow directly for any 
radius: Pg(z) = pg(Rj^)z“-^, Pg(z) = Xg(z) Pq(Rj^)z“‘^; v^, Fg, and Ug 

follow from equation (1). The quantity ~ Xg(2.53-(u) can then 

be used in the two equations (3) to find a and $. Equation (13) is 
used to find n and v follows from the identity ot-0 = \M-ri. Your 
algebra can be checked by comparing the results with those found from 
equation (7). This procedure supplies the postshock conditions and 
their derivatives at any value of z. The complete density structure is 
provided by equation (20). In order to find the remainder of the 
structure or the time evolution of an individual parcel, one performs 
numerical integrations as in Section II. 

Suppose one wishes to learn the past history of a particular gas 
parcel now located at r^ =» R^/Rg. Equation (20) provides its current 
mass fraction pj and at any other epoch, y = y£(Zp/z)^”‘^. Thus the 
parcel was first encountered by the shock when z = z^^ = Zp( 
when y = y^ =* !• Conditions in the parcel at that moment are found as 
before since it is then the postshock gas. From these conditions one 
evaluates = [xg(zj) and = [l.l/1.2]^^^ K^. 

From that point, one steps forward in z, following the time with 
equation (19) and the evolution of K and C with equations (21). At 
each new shock position z, one must calculate y for the parcel, invert 
equation (20) to find r and then x (using the new postshock values of 
Xg and a), find n = xnQ(R£)z~‘^, calculate T and g = T^/T from C and K, 
calculate (if desired) p from n and T, calculate any desired 
ionization and recombination rates and derivatives of ionic 
concentrations, and finally calculate new derivatives of K and C from 
equations (22), (23), and (24) using the contemporary values of K, C, 
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g, X, y, x' and y for the parcel and of z, x^, v at the shock. The 
procedure appears somewhat messy but it adds very little to the number 
of integrals which have to be performed anyway to follow the evolution 
of the ion concentrations. 


IV. Discussion 

The present work shows that including the effects of saturated 
thermal conduction depresses the shock compression factor x^ only 
modestly from the value of 4 appropriate to adiabatic calculations. 
We also find that the solution undergoes a sharp transition from the 
early time asymptote where thermal conduction dominates to the late 
asymptote where the structure is essentially adiabatic and the 
electron and ion temperatures are nearly equal over a large fraction 
of the remnant. The radius at which this occurs is quite similar to 
that given by Cowie (1977) for the equilibration of the two 
temperatures. 

An error in the captions of Figures 2 and 3 of Paper I should be 
pointed out: actually the ambient density should be a Rg^ and PgRg^ 
respectively. 

This work was supported in part by NASA grant NGL 50-002-044 at 
the University of Wisconsin-Madison. 
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Table 1 

Variables, Derivatives, and Normalizations 


Symbol Postshock Derivative 

Variable f ( 3£nf/ 3JlnR)g Normalized Form 


radius 

R 


r = 

R/R^ 

shock radius 

R 

s 


z = 

Rg/Rl 

time 

t 


T = 

t/ti 

density 

P 

a= ( 31.np/3lnR)g 

X = 

p/Pq(Rs) 

nuclear number 
density 

n = 

p/(1.4 mjj/1.1) 



pressure 

P 

0 = ( 31np/ 3l.nR)g 

y = 

p/ps 

average temperature 

T 

(0-a) = (3ilnT/3ilnR)^ 

s 



electron temperature 

^e 

n = -(31.nTg/3AnR)g 



mass velocity 

u 

0 = ( 3tou/ 3£nR) 

5 



shock velocity 





thermal conduction 
flux 

F 




div F 

1 
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Table 2 

Solution Properties 


g(u))^ 


(1) 

z® 

max 


Early 

Late 

0 

.858 

3.2383 

0.643 

0.655 

-2 

.970 

2.9121 

0.578 

0.581 

-4 

.998 

2.7721 

0.511 

0.549 


^the value of z where x* Is maximum 


^e(w) = Eq/ [ 3VgPg/2 ]; from Paper 1 
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Table 3 

Analytical Fit Coefficients 


(0 


Cl 

C 2 

C3 

dl 

^2 

83 

0 

3.2383 

.37903 

.069789 

.014860 

.32173 

.059507 

.012030 

-2 

2.9121 

.50027 

.11527 

.040463 

.38703 

.091284 

.029458 

-4 

2.7721 

.54857 

.13277 

.054586 

.40554 

.10147 

.037829 



-24- 


References 

Cowle, L. L. and McKee, C. F. 1977, Ap. J. , 2ll, 135. 

Cowie, L. L. 1977, Ap . J . , 215, 226. 

Cox, D. P. and Anderson, P. R. 1982, Ap . J . , 253, 268 (CA) . 

Cox, D. P. and Edgar, R. J. 1983, Ap ♦ J . , 265, 443 (Paper I). 

Cox, D. P. and Franco, J. 1981, Ap. J. , 251, 687 (CF). 

Kahn, F. D. 1975, in Proc. Int. Cosmic Ray Conf. , (Munich), 11, 3566. 


-25- 


Figure Captions 

Fig. 1 - Values of edge parameters and derivatives plotted as 

functions of normalized shock radius z=R /R. , for ambient densities 

s ^ 

Po“^r 

a) The compression factor h Pg/pQ 

b) X* = 3ilnXg/3nRg 

c) The logarithmic postshock density derivative 

a = (3inp/3inR)g 

d) The logarithmic postshock pressure derivative 

3 = ( 3Jinp/ 3inR) ^ 

e) The logarithmic postshock electron temperature derivative 
n = -(3£nTg/3JlnR)g 

f) The logarithmic derivative of the electron-to-average 
temperature ratio v = [9Jln(Tg/T)/ 3JlnR]g 

Fig. 2 - Normalized age of the remnant T = t/t^^ as a function of 
noirmalized shock radius z = Rg/R]^ 

a) 0 ) = 0, a uniform ambient density 

b) Pq“R2 

c) Pq“R*^ 

Fig. 3 - Snapshots of the internal remnant structure for 
Pq = constant 

a) Pressure in units of Pg(Rj^) 

b) Average temperature (solid lines) and electron 


temperature in units of Tg(Rj^) 
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c) Density in units of Pq(R^)- Also plotted are ambient 

density (lower line) and post-shock density (upper line). 

9 ■' 

Fig. 4 - Same as Figure 3 but with p “ R„. 

o s 

Fig. 5 - Same as Figure 3 but with Pq “ Rg* 
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